Transmission through Quantum Dots: Focus on Phase Lapses 
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Measurements of the transmission phase in transport through a quantum dot embedded in an 
Aharonov-Bohm interferometer show systematic sequences of phase lapses separated by Coulomb 
peaks. Using a two-level quantum dot as an example we show that this phenomenon can be ac- 
counted for by the combined effect of asymmetric dot-lead coupling and interaction-induced "pop- 
ulation switching" of the levels, rendering this behavior generic. In addition, we use the notion of 
spectral shift function to analyze the relationship between transmission phase lapses and the Friedel 
sum rule. 

PACS numbers: PACS numbers: 73.21.La, 73.63.Kv, 73.23.Hk, 03.65.Vf 
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In a series of experiments by the Weizmann group, 
the transmission phase, O tr , characterizing transport 
through a quantum dot (QD) has been systematically 
studiedi*2»2, embedding the QD in an Aharonov-Bohm 
interferometer^. Arguably the most intriguing finding 
of these experiments has been the correlated behavior 
of Q t r as function of the leads' chemical potential fi (or 
the gate voltage): it appears to undergo a lapse (phase 
lapse, PL), seemingly of — n, between any two consecu- 
tive Coulomb peaks. It is clear that this effect cannot 
be explained within a single-particle framework™. More- 
over, in spite of a substantial body of theoretical work 
(see, e.g., Refs. [Hi), some of which gained important 
insight on the underlying physics, no clear cut theory- 
experiment connection has been established as yet. 

In the present article we revisit this problem. We do 
this by studying a (spinless) two- level QD, attached to 
two leads. We account for the difference in the couplings 
of level 1 and level 2 to the leads ("1 — 2 asymmetry" ) 
and, for the first time, probe the effect of the (generi- 
cally expected) asymmetric coupling to the left and the 
right leads ("L-R asymmetry"). We find unexpectedly 
that these two asymmetries give rise to a qualitatively 
new behavior of tr (/i), and render the appearance of 
PL between consecutive Coulomb peaks generic. This 
conclusion is in line with recent renormalisation group 
results for a QD with degenerate levels^. 

Throughout the discussion of transmission PLs in the 
literature, much attention was paid to the Friedel sum 
rule, which, in one dimension, relates the transmission 
phase to the change of carrier population in the sys- 
tem (see, e.g., Refs. [ITIfT^ ). Since the latter varies 
monotonously with the chemical potential (or gate volt- 
age) , one may perceive a contradiction between this sum 
rule and the occurrence of PLs. We revisit this issue in 
Appendix [Al and show, in particular, that the correct for- 
mulation of Friedel sum rule in one dimension allows for 
transmission phase lapses. 

The minimal model for studying the phase lapse mech- 
anism includes a two-level QD, 



Here, the operators di with i = 1,2 annihilate elec- 
trons on the two dot sites (with bare energies 25- , 



E { 2 0) > E\ v >). The QD is coupled to the two leads by 
the tunnelling term 
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The operators Cj (with half-integer j) are defined on the 
tight-binding sites of the left and right lead (cf. Fig. [1]). 

We begin with summarizing the results of Refll3l (see 
also Refs. [l2lfl3ll5llT6h in the case when no charging in- 
teraction is present, U = 0, and the value of G tr is readily 
calculated (even for a larger number of dot levels). The 
two transmission peaks then take place near /i = , 
each corresponds to a smooth increase of 8t r (^) by tt 
within a chemical potential range proportional to aj^+a 2 ^ 
for the first dot level, b\ + b\ for the second one. If 
the relative coupling sign, a = sign(a£OR&L&.R), equals 
+1 (same-sign case), a discontinuous PL of A0 tr = — tt 
(transmission zero) arises in the energy interval between 
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H QD = (E[ 0) - fJl )d\d 1 + (E [ 
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the two transmission peaks, E^ < ji < E^ . 
this would be in qualitative agreement with the mea- 
surements, experimentally there is no way to control the 
coupling signs. Indeed, for the relevant case of a random 
(chaotic) QD, one expects close to 50 % of the adjacent 
pairs of dot levels to have a = — 1 (opposite-sign case), 
when no phase lapse occurs between the two correspond- 
ing level crossings. These observations^ (and hence the 
difficulty in accounting for the experimentally observed 
correlations in Qtr) persist even when interaction is ac- 
counted for (but when |ol| = |or| = \bh\ = \bn\ was 
assumed) . 

Following the original idea of Ref. [l7|, the effects of 
"population switching" due to a charging interaction U 
in discrete spectrum QDs [Eq. {I])] were addressed both 
theoreticall y 18 ' 19 ! 20 ' 21 ! 22 and experimentally 2 ^,. If one of 
the dot levels is characterized by a stronger coupling to 
the leads and U is sufficiently large, the two level oc- 
(d\di) show non-monotonic dependence 
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on /i. A rapid "population switching" 20 i 21 (which may 
be accompanied by the switching of positions of the two 
mean-field energy levels, £1,2), takes place. The avail- 
able results, however, remain incomplete in that (i) the 
behavior of rii near switching (abrupt vs. continuous for 
different values of parameters) was not investigated, (ii) 
only the case of \cll\ = \or\ and = |&r| was con- 
sidered, omitting the important effects of coupling asym- 
metry (see, however, Ref. [l(J), and (iii) the relationship 
between population switching and PLs was not addressed 
fully and correctly. The present article is aimed, in part, 
at clarifying these issues. 

We find that at sufficiently large U, including the 
dot-lead coupling asymmetry largely alleviates the "sign 
problem" as outlined above, giving rise to a phase lapse 
of A0 tr = — 7r between the two Coulomb peaks for the 
overwhelming part of the phase diagram at both a = 1 
and o~ = — 1. This is a result of an effective renormaliza- 
tion of the coupling sign, er = — 1 to a = 1, due to the 
interaction. As some asymmetry of individual level cou- 
pling is generally expected in experimental realizations of 
QDs, this novel phase-lapse mechanism appears relevant 
for understanding the experimental data. Furthermore, 
we consider the implications of interaction-induced "pop- 
ulation switching" 17 ! 18 for the transmission phase. We 
show that, under certain conditions ("abrupt" switch- 
ing), this leads to a modification of phase-lapse value 
(|AO tr | < 7r). Once fluctuations (omitted in the present 
mean-field treatment) are taken into account, this result 
may translate into a more complex behavior in the vicin- 
ity of the phase lapse. 

The analysis of the full four-dimensional space of all 
values of ol.r and proves too cumbersome and per- 
haps redundant. Rather, we find it expedient to investi- 
gate a suitable 3D subspace, which is defined by a con- 



straint, b R 



i R . Then there exists a unitary 



transformation of the two dot operators, d\^ 2 — * ^1,2, 
changing the coefficients in Eq. ([2]) in such a way that 
o-l = o,r = a, 6l = ahji = b with a = — 1 (the a = 1 case 
corresponds to the same-sign symmetric original cou- 
pling: cll = &l = bn). The transformation also 
affects the form of the first two terms on the r. h. s. of 
Eq. ([T]), which now read 

(^ 0) - v)d\di + (Ei 0) - M )4<i~2 - ^-(d\d 2 + ^di). (3) 

The coefficients e[°1 and wq can be formally thought of 
as the bare "site energies" and "intra-dot hopping" of 
a QD depicted in Fig. [TJ and are related to the level 
energies [cf. Eq. ©] by 2E[ 0) 2 = (E^ + E^ ) =p [(e[ 0) - 

I?2 ) 2 + w o 2 ] ■ Our analysis will be carried out in terms 
of this new QD with a = ±1. For the a = — 1 case, u>o 
is actually a measure of (left-right) asymmetry in the 

coupling of the original QD levels, E\ °l , to the two leads. 

Our calculation consists of the following steps: (i) mean 
field decoupling of the interaction term in[TJ (ii) obtain- 
ing an effective single particle Hamiltonian in terms of 
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FIG. 1: The model system, composed of a wire (chain) and a 
two-level dot, Eqs. and (J3J). 



the averages (d\dj); (iii) expressing O tr in terms of the 
parameters of that Hamiltonian; (iv) expressing {d\dj) in 
terms of Qtr employing the Lifshits-Krein trace formal- 
ism; (v) solving self-consistently the resultant equations 
for {d\dj) ; (vi) obtaining explicit results for Qtr- 
(i) Mean field decoupling reads 



d\dld 2 di -> d\di{d\d 2 ) + dld 2 (d[di) - 



i«i 



)(d\d 2 ) 



-dld 2 (d 2 di) - 4^1(^1^) + \{d\d 2 )\ 2 . (4) 

We verified that the results of our mean-field scheme are 
independent on the basis (of the two dot states) in which 
the decoupling is carried out. In the case of asymmet- 
ric coupling, it is important^ to keep the off-diagonal 
(" excitonic" ) average values in the above expression, e.g. 
{d\d\). Owing to a cancellation between virtual hopping 
paths between the two QD sites, these averages vanish in 
the er = — 1 symmetric case of wq = (corresponding to 
«l = o,r, bi = — bn) 20 i 21 . However, this does not occur 
generally, nor indeed in the same-sign symmetric case, 
leading to difficulties noted in Ref. 

(ii) Substituting Eqs. ([SHU into |T]) is tantamount 
to mapping of the original model onto an effective non- 
interacting model with the Hamiltonian given by Eq. (|2|) 
and the mean-field dot term, 

Hf F = {E x - ii)d\dx + {E 2 - n)d\d 2 - ^(d\d 2 + 4^). 

(5) 

The self-consistency conditions take the form of three 
coupled mean-field equations, 



E 1 = El 



(0) 



U (did 



2"2/, 



E2 — E, 



(0) 



w = w + 2U(d\d 2 ). 



U{d\d 1 ), (6) 
(7) 



(iii) For the effective single-particle model §5$ one can 
readily compute the transmission phase, Ot, (e). In the 
a = — 1 case, it is given by 



fT 2 2+ n , b 2 (E 1 -e) + a 2 (E 2 -e) + 2e± 

V t 2 - e z tan9 tr = e H Tl 

[E 1 -e){E 2 -e)- 1 ^ ^ 



i w2 ~ s 



b- 

b 2 

(where 2t is the width of conduction band in the leads) 
and suffers a lapse of — ir at that value of e for which the 
transmission vanishes, i. e., e — Z, 



E x b 2 



b 2 



(9) 
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(iv) The quantum mechanical average values in Eqs. 
([5H7]) are given by derivatives 

(d\ 2 Ji.2> = dQ M F/dEi,2 , (d\d 2 ) = -dVL MF /dw 

( 10 ) 

of the thermodynamic potential of the effective system. 
The latter is evaluated exactly with the help of the 
Lifshits-Krein trace formula^ 5 -, 



MF 



ae)de. 



(11) 



Here, O is the combined potential of a disconnected sys- 
tem comprising a dot [Eq. (J3J)] and a wire, 



(12) 



j j 
The spectral shift function £ is defined by its relationship, 

d^/de = -5v(e) , (13) 

to the change of the total density of states of this system 
due to a local perturbation, 



V = V T + t - (c\ /2 c_ 1/2 + c T _ 1/2 c 1/2 ) 



(14) 



In Eq. ([14]). the second term on the r.h. s. corresponds 
to cutting the link between sites j = —1/2 and j = 1/2 
of the wire. The two resulting leads are coupled to the 
QD by Vt [Eq. ||3J|]. Since for a wire of a finite length 
2L, £ is related to the shifts of (discrete) energy levels 
under the effect of V, it is easy-ii to express £ in terms 
of ®tr, viz. 



-6 tr /7r + m(e) 



(15) 



(see Appendix |XJ). The integer- valued function m should 
be chosen to satisfy the requirement 25 for £ (e) to vanish 
continuously with decreasing strength of the perturba- 
tion (e.g., XV with A — * 0). We find that the value of m 
changes by +1 at e = E lj2 [eigenvalues of TC^ F , Eq.©], 
and by —1 at the transmission zero. For a — — 1 we 
obtain (in the units where t = 1): 



(did!) = 



(d\d 2 ) 



VI - e 2 X(e)de x (16) 

-l 

a 2 b\l-e 2 )+a 2 (E 2 -e + b 2 ef+ b ^ 
Vl- e 2 X(e)dex (17) 



a% 2 (l-e 2 )+b 2 (E x -e + a 2 e) 2 + 



(d\d 2 



V 7 ! - e 2 X(e)de x 



(18) 



a 2 (E 2 - e) + b 2 (Ei - e) + 2a 2 & 2 e 



TrX(e) 



(1-e 2 ) a\E 2 -e)+b 2 (E l -e) + 2a 2 b 2 e 



-i 2 



(£i -e + a 2 e)(£:2 - e + 6 2 e) - (1 - e 2 )a 2 6 2 - -w z 

Similar expressions are obtained also for the a = 1 case 
(see Appendix IB"]) . 

We now sofoe equations ((Mil) numerically (v), and sub- 
stitute the resulting values of Ei, 2 and w into the expres- 
sion for tr to get the transmission phase (vi). 

The results are summarized in the phase diagram, Fig. 
[2] The left-hand part corresponds to a — 1, whereas 
the a = — 1 case (when the results do not depend on 
the sign of wq) is shown on the right. The bold line 
marks the boundary between continuous (phases 1-2) 
and discontinuous (see below) regimes of dependence of 
the effective QD parameters on /i. Within each regime, 
different phases are identified according to the magni- 
tude and location of PL(s) with respect to the trans- 
mission peaks (Fig. [2] table). It should be noted that 
in the cr = — 1 case the latter are given by /ii : 2 = 
(E t + E 2 )/2 T - E 2 ) 2 + w 2 + \a 2 b 2 lt 2 \ x l\ and 

are slightly shifted with respect to mean field dot levels, 
£i, 2- In the table, we denote transmission peaks by E\ t2 
irrespective of the sign of a in order to keep the notation 
uniform. Typical dependence of Qtr on /x for each phase 
is shown in Fig. [3J 

In the continuous-evolution part, phase 1 (phase 2), 
which occupies a large (small) area of the phase space, 
corresponds to the case when the phase lapse of — tt, as- 
sociated with the transmission zero, lies within (outside) 
the interval of values of /x between the two transmission 
peaks. It should be noted that the right-hand, a = — 1, 
side is expected to be representative of both opposite- 
and same-sign cases (a = ±1), provided that the left- 
right asymmetry is sufficiently strong (large Wq). This 
is illustrated by the thin solid line, above (below) which 

coupling signs for the two bare dot levels E^l become 
the same, a = 1 (opposite, a — —1). Once the in- 
teraction effects are taken into account, one sees that 
phase 1 extends also far below this line, which is indica- 
tive of the effective change of the coupling sign [due in 
turn to the interaction-induced enhancement of w; at 
w ^> E 2 — Ei, the coupling of the two mean-field dot 
levels, (di ± d 2 )/V2, to the leads is same-sign]. 

The discontinuous behavior is associated with the pres- 
ence of multiple solutions of the mean field equations JS]- 
[7]) within a range of values of /z, which is illustrated by 
a "fold" (bold solid and dashed lines) on the schematic 
Ei(ijl) — [i plots in Fig. [3] We find that if a system for- 
mally is allowed to follow such a multiple- valued solution 
from left to right, the value of Qtr increases, and also 
suffers a PL of — n at some point (marked by a circle). 
In reality, thermodynamics dictates that the full thermo- 
dynamic potential = Qmf — U (d\d~i) (d\d 2 ) + U (d\d 2 ) 2 
[cf. Eq. (jTTj) ] should be minimized to identify the stable 
solution, resulting in a "jump" (vertical line), which in 
turn is associated with a positive increase of fr by a frac- 
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- transmission peaks; Z; - phase lapse locations) 



FIG. 2: (color) The "phase diagram" of a two-level QD with Sl = an — a, bh = abR — b. The parameters are U = O.lt, 
E{ 0) = 0, Ef ] = 0.004t, and Va 2 + b 2 = 0.125*. The axes represent the 1-2 level asymmetry, a = (\a\ - \b\)/^a 2 +b 2 , and the 
dimensionless intra-dot hopping, n = wo /[(E[ 0) - Ef y ) 2 + w' 2 ] 1/2 . Properties of different phases are summarized in the table. 
At U < 0.04t, the border of discontinuous-evolution region (bold line) does not meet the boundary between phases 1 and 2. 
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FIG. 3: (color) Typical behavior of 6t r (A0 in different phases 
(top left; plots shifted for convenience). Relative positions of 
transmission peaks (boxes; also in the main panel) and the 
— 7r-PLs (circles) in phases 3-7 are clarified by the schematic 
plots of Ei — n around the multiple-solution region (absent for 
phases 1-2). Solid (dashed) lines correspond to stable (unsta- 
ble) solutions. The abrupt "switching" of solutions (vertical 
solid line) may either renormalize the PL (when the — 7r-lapse 
lies in the unstable region) or result in a positive jump of Qtr- 



tion of 7r, giving rise to a second "PL" (phases 4,6), and 
with the population switchin g 17 i 18 of the dot levels. If the 
transmission zero lies in the thermodynamically unstable 
part of the solution (bold dashed line; phases 3,5,7), the 
PL of — 7r should be added to this positive increase of 
Qtr, giving rise to a single "renormalized" PL. Finally, 
one of the transmission peaks may be located within the 
unstable region (phases 5,6) with a result that the plot 
of &tr(p) does not have a corresponding inflection point, 
which is replaced by a PL. 

It follows that at least within the mean-field frame- 
work discontinuous population switching is always asso- 
ciated with the presence of multiple solutions and hence 
with "renormalized" PLs (or alternatively with addi- 
tional "PLs" characterized by an increase of phase by a 
fraction of 7r). This conclusion is clearly at variance with 
the suggestion of Ref. [T7| that the discontinuous switch- 
ing between multiple solutions gives rise to the PLs of 
7r as observed experimentally. We note that while the 
behavior of transmission phase in this regime should be 
investigated beyond the mean field, the main point of our 
paper is that there is another mechanism which gives rise 
to a PL of 7T without a discontinuous population switch- 
ing (phase 1). Since this latter scenario does not involve 
instabilities of any kind, it can be expected to remain 
robust with respect to fluctuations (not included in the 
present treatment). 

In summary, we have presented here a generic mech- 
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anism for the appearance of phase lapses between 
Coulomb blockade peaks. These PLs may be renormal- 
ized by a discontinuous "population switching". Exper- 
imentally it would be interesting to correlate the latter 
with the former by simultaneously measuring dot occu- 
pancy (employing a quantum point contact), and trans- 
mission phase. Theoretically, going beyond a mean field 
analysis is needed to determine the importance of quan- 
tum fluctuations. 
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APPENDIX A: SPECTRAL SHIFT FUNCTION, 
TRANSMISSION PHASE, AND FRIEDEL SUM 
RULE 

For the case at hand, the use of the standard formula 25 
for the spectral shift function £, 



£(e) = -iArgDct (l — — 1 - 

7T £ — lO — ft,. 



V 



(Al) 

proves rather cumbersome. Instead, we will use the un- 
derlying notion of spectral shifts^ 5 - in order to derive 
the generic relation (|15[) between £ and the transmis- 
sion phase. This derivation also allows for an important 
insight concerning the Friedel sum rule. 

We consider a system similar to that shown in Fig. [TJ 
with the QD between the sites —1/2 and 1/2 replaced 
by an arbitrary point scatterer. The latter is charac- 
terized by an S'-matrix whose elements have a smooth 
dependence on the particle energy 2 ^. While the bound- 
ary conditions cannot affect the value of £ in the limit 
when the length of the wire, 2L, is large, the treatment is 
simpler when periodic boundary conditions are assumed. 
The spectrum of the wire in the absence of the scatterer, 
which we refer to as unperturbed, is then given by 



e(kj) = — cos kj 



k - H 



j=0,l,2,...,L (A2) 



[cf. Eq. (fl"2|) where we assumed t = 1]. The wave func- 
tions are proportional to exp (±ikjx) and, for j ^ 0, L, 
the corresponding energy levels are doubly degenerate. 
Since we are ultimately interested in the L — > oo limit, it 
is assumed that the inter-level spacing in the wire consti- 
tutes the smallest energy scale in the problem. The levels 
are shifted, and the degeneracy is lifted, in the presence 
of the scatterer, when the wave function is generally given 
by 



A 2 e ikx + B 2 e 



-\kx 

) 

-\kx 



x < 0, 
x > 0. 



(A3) 



The linear relationship between coefficients on the right 
and on the left of the scatterer reads (assuming time- 
reversal symmetry 



A 2 = aAi + /3Bi , B 2 = (3*Ai + a*B x 



(A4) 



with \a\ 2 — |/3| 2 = 1. Relation of the quantities a and 
(3 to the ^-matrix is given by, e.g., setting Ai = (in- 
coming particle from the right), hence (right-right) reflec- 
tion amplitude, r rr — 13/ a* and transmission amplitude, 
t tr = 1 /a*. 

Now the periodic boundary conditions dictate that the 
allowed momentum values shift, 



(12) Aj^TT 



(A5) 



Substituting Eqs. (|A3HA5[) into the condition %p{—2L 
0) = ip(+0), we find for 1 = 1,2 

Aie -27riA< ! > = A 2 = aA 1 +(3B 1 , 
Bie 27riAW = B 2 =l3*A 1 +a*B l . 



This yields the equation for A^- 1 ' 2 ^ (cf. Ref. [Til ): 



ae 4 - iA(i> - 2e 2wiA(0 + a* = 0, 



(A6) 



or equivalently 



^iA<') _ e 2-iA< 1 )^ ^iA<" _ e 2.iAf) 



2 4.iA<0 _ 2 e 

a 



27riA 



(0 a 
a 



yielding 



^ifA^+Af) = or 



(A7) 



(1 2) 

In the limit L — > oo, the quantities A^- ' become func- 
tions of energy and, writing also a = exp(i0 tr )/|i tr | with 
Qtr the transmission phase, we find 

A[e(fc J )] = Af+Af = —QtrWkM+mlefa)] . (A8) 

Let us now discuss the quantities appearing in Eq. 
(IA8[) . (i) Qtr is the transmission phase. In the pres- 
ence of localized states within the scatterer (dot levels 
Ef), Qtr increases by n as the energy of interest [e(kj) in 
our notation, or more physically, the chemical potential] 
spans a resonance, (ii) m(e) is an integer which we will 
now choose in such a way that A coincides with the Lif- 
shits - Krein spectral shift function, £(e) [see Eq. (fT3|) ]. 
m(e) then changes by +1 with increasing energy at ev- 
ery Ef; in addition, it changes by —1 at the points where 
transmission vanishes (transmission PLs). We thus arrive 



at Eq. |H]). (iii) Ay'ir/L and A) >ir/L are the shifts in 
the allowed values of momentum [cf Eq. (|A5|> ]. 
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There is no bound state corresponding to a PL, imply- 
ing that £(e) should be continuous at that point (trans- 
mission zero). The choice of m(e) discussed above [along 
with Eq. (|A8[) ]. ensures that A indeed may be synony- 
mous with £ (see below). 

In order to use the calculated value of spectral shift 
function for the total energy evaluation via the trace 
formula, Eq. (fTTj) . one needs to know the overall ad- 
ditive constant in £(e). In the regime of interest to us, 
no bound state is formed below the band bottom (at 
e < —1). From the viewpoint of the lowest-energy elec- 
tron states (e — > —1 + 0), the scatterer then acts as an 
impenetrable potential barrier (and not as a potential 
well), and the constant is fixed by a readily derivable 
condition, £(e — > —1 + 0) = 1/2, valid for any barrier 
with no bound state formed below its bottom. 

The quantity £(e) remains a smooth function of energy 
e away from band edges e — ±1 and the dot levels Ef. 
As mentioned in the text, the spectral shift function is 
related to the perturbation-induced change in the density 
of states. For the unperturbed system, the latter can be 
defined in the L — > oo limit only as 



KM] = I™ 



eOj+i) - e(^) 



oc L . 



(A9) 



where the factor of 2 reflects the double degeneracy of 
energy levels. In the presence of the scatterer we obtain, 
with the help of Eq. JA5|, 



where the first term on the r. h. s. is the band filling 
of an unperturbed wire. By re-writing this in terms of 
transmission phase 6 tr (e) [cf. Eq. (|A8[) ], we g e * the 
Friedel sum rule in the form 



(A12) 



With increasing /i, the integer to changes by +1 at trans- 
mission zeroes, /i = Zi. We note that the sum of the two 
last terms on the r. h. s. of Eq. (|A12[) remains contin- 
uous at fi = Zi, emphasizing that the Friedel sum rule 
does not account for the transmission phase lapses. This 
is because the underlying spectral characteristic, £(e) [cf. 
Eq. IjAlip ] remains smooth at /i = Zi and in general 
does not depend on \t tr \- 



APPENDIX B: MEAN-FIELD EQUATIONS IN 
THE 5 = 1 CASE 



In the case of same-sign symmetric coupling of the QD 
to the leads, a = 1, Eqs. ([HM]) and ([161 [PS]) are replaced 
with 



\Jt 2 — e 2 tanOt,- = e + 



& 2 (£i - e) + a 2 (E 2 - e) + abw 



(E 1 -e)(E 2 -e)-~w 



1 2 ' 



<m = 

lim 



[e(k j+1 ) - eikj)} {1 + ±£[e(fc J+1 )] - 
= u [e(kj)] + Su[e(kj)] , 



where 



Sv(e) = -d£(e)/de. 



(AlO) 



Here, we used the obvious fact that the centre of grav- 
ity of the two perturbed levels formed out of a doubly 
degenerate unperturbed level e(fcj) is given by e(kj) + 
\[e{kj+i) — e(kj)]A[e(kj)] (substituting A with £ for our 
choice of to). We note that, as expected on physical 
grounds, the quantity 5u is not extensive, i. e., it is not 
proportional to the length of the wire (in contrast to vq). 
For the specified choice of m(e) in Eq. (|A8[) . Eq. (|A10[) 
yields also a delta-functional contribution to 5v of the 
form — J2i — Ef). This corresponds to merging of the 
discrete dot levels into continuum and shows that 5v is 
the difference in the density of states between the wire 
with the scatterer and a disconnected system comprised 
of an unperturbed wire alongside an isolated scatterer. 

Integrating Eq. (|A10[) . we get the expression for the 
total particle number, 



and 



{d\dx) = 



(d\di) = 
(d\d 2 ) = 



Z = 



E 2 a 2 + Eib 2 + abw 
a 2 + b 2 



1 2 



VT 



a{E 2 — e) + -bw 



b(E\ — e) + 



Y(e)de, 



Y(e)de, 



a/1 - e 2 Y(e)de x 



a(E 2 — e) + -bw 



b{E\ — e) + 



respectively. Here, 
1 



= (l-O 



( El -e)(E 2 -e)--w 2 



+ <j a 2 (E 2 - e) + b\E 1 - e) + abw + 
(S 1 -e)(£; 2 -e)- iU ; 2 
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